On neoclassical impurity transport in stellarator geometry

The impurity dynamics in stellarators has become an issue of moderate concern due to the inherent tendency of the impurities to accumulate in the core when the neoclassical ambipolar radial electric field points radially inwards (ion root regime). This accumulation can lead to collapse of the plasma due to radiative losses, and thus limit high performance plasma discharges in non-axisymmetric devices.\\ A quantitative description of the neoclassical impurity transport is complicated by the breakdown of the assumption of small $\mathbf{E}\times \mathbf{B}$ drift and trapping due to the electrostatic potential variation on a flux surface $\tilde{\Phi}$ compared to those due to the magnetic field gradient. The present work examines the impact of this potential variation on neoclassical impurity transport in the Large Helical Device (LHD) stellarator. It shows that the neoclassical impurity transport can be strongly affected by $\tilde{\Phi}$. The central numerical tool used is the $\delta f$ particle in cell (PIC) Monte Carlo code EUTERPE. The $\tilde{\Phi}$ used in the calculations is provided by the neoclassical code GSRAKE. The possibility of obtaining a more general $\tilde{\Phi}$ self-consistently with EUTERPE is also addressed and a preliminary calculation is presented.Comment: 11 pages, 15 figures, presented at Joint Varenna-Lausanne International Workshop on Theory of Fusion Plasmas, 2012. Accepted for publication to Plasma Phys. and Control. Fusio

Comparison of particle trajectories and collision operators for collisional transport in nonaxisymmetric plasmas

In this work, we examine the validity of several common simplifying assumptions used in numerical neoclassical calculations for nonaxisymmetric plasmas, both by using a new continuum drift-kinetic code and by considering analytic properties of the kinetic equation. First, neoclassical phenomena are computed for the LHD and W7-X stellarators using several versions of the drift-kinetic equation, including the commonly used incompressible-ExB-drift approximation and two other variants, corresponding to different effective particle trajectories. It is found that for electric fields below roughly one third of the resonant value, the different formulations give nearly identical results, demonstrating the incompressible ExB-drift approximation is quite accurate in this regime. However, near the electric field resonance, the models yield substantially different results. We also compare results for various collision operators, including the full linearized Fokker-Planck operator. At low collisionality, the radial transport driven by radial gradients is nearly identical for the different operators, while in other cases it is found to be important that collisions conserve momentum

Properties of a new quasi-axisymmetric configuration

A novel, compact, quasi-axisymmetric configuration is presented which exhibits low fast-particle losses and is stable to ideal MHD instabilities. The design has fast-particle loss rates below 8\% for flux surfaces within the half-radius, and is shown to have an MHD-stability limit of a normalised pressure of $\langle\beta\rangle=3\%$ where $\langle\beta\rangle$ is volume averaged. The flux surfaces at various plasma betas and currents as calculated using the SPEC equilibrium code are presented. Neoclassical transport coefficients are shown to be similar to an equivalent tokamak, with a distinct banana regime at half-radius. An initial coil design study is presented to assess the feasibility of this configuration as a fusion-relevant experiment

Simulating Gyrokinetic Microinstabilities in Stellarator Geometry with GS2

The nonlinear gyrokinetic code GS2 has been extended to treat non-axisymmetric stellarator geometry. Electromagnetic perturbations and multiple trapped particle regions are allowed. Here, linear, collisionless, electrostatic simulations of the quasi-axisymmetric, three-field period National Compact Stellarator Experiment (NCSX) design QAS3-C82 have been successfully benchmarked against the eigenvalue code FULL. Quantitatively, the linear stability calculations of GS2 and FULL agree to within ~10%.Comment: Submitted to Physics of Plasmas. 9 pages, 14 figure

Translating evidence into practice : ACOs’ use of care plans for patients with complex health needs

Background Care plans are an evidence-based strategy, encouraged by the Centers for Medicare and Medicaid Services, and are used to manage the care of patients with complex health needs that have been shown to lead to lower hospital costs and improved patient outcomes. Providers participating in payment reform, such as accountable care organizations, may be more likely to adopt care plans to manage complex patients. Objective To understand how Medicare accountable care organizations (ACOs) use care plans to manage patients with complex clinical needs. Design A qualitative study using semi-structured interviews with Medicare ACOs. Participants Thirty-nine interviews were conducted across 18 Medicare ACOs with executive-level leaders and associated clinical and managerial staff. Approach Development, structure, use, and management of care plans for complex patients at Medicare ACOs. Key Results Most (11) of the interviewed ACOs reported using care plans to manage care of complex patients. All care plans include information about patient history, current medical needs, and future care plans. Beyond the core elements, care plans included elements based on the ACO’s planned use and level of staff and patient engagement with care planning. Most care plans were developed and maintained by care management (not clinical) staff. Conclusions ACOs are using care plans for patients with complex needs, but their use of care plans does not always meet the best practices. In many cases, ACO usage of care plans does not align with prescribed best practices: ACOs are adapting use of care plans to better fit the needs of patients and providers

Study of the neoclassical radial electric field of the TJ-II flexible heliac

Calculations of the monoenergetic radial diffusion coefficients are presented for several configurations of the TJ-II stellarator usually explored in operation. The neoclassical radial fluxes and the ambipolar electric field for the standard configuration are then studied for three different collisionality regimes, obtaining precise results in all cases

The monoenergetic approximation in stellarator neoclassical calculations

In the standard "monoenergetic" approach to numerical calculation of stellarator neoclassical transport, to expedite computation, ad-hoc changes are made to the kinetic equation so speed enters only as a parameter. Here we examine the validity of this approach by considering the effective particle trajectories in a model magnetic field. We find monoenergetic codes systematically under-predict the true trapped particle fraction, with the error in the trapped ion fraction being of order unity when the electric field is large, suggesting some results of these codes may be unreliable in this regime. This inaccuracy is independent of any errors introduced by approximation of the collision operator.Comment: 8 pages, 2 figures, submitted to Plasma Phys. Controlled Fusio